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Abstract We study the global regularity, for all time and all initial data in 
H 1 / 2 , of a recently introduced decimated version of the incompressible 3D 
Navier-Stokes (dNS) equations. The model is based on a projection of the 
dynamical evolution of Navier-Stokes (NS) equations into the subspace where 
hclicity (the L 2 — scalar product of velocity and vorticity) is sign-definite. The 
presence of a second (beside energy) sign-definite inviscid conserved quadratic 
quantity, which is equivalent to the H 1 / 2 — Sobolev norm, allows us to demon- 
strate global existence and uniqueness, of space-periodic solutions, together 
with continuity with respect to the initial conditions, for this decimated 3D 
model. This is achieved thanks to the establishment of two new estimates, 
for this 3D model, which show that the H 1 / 2 and the time average of the 
square of the H 3 / 2 norms of the velocity field remain finite. Such two addi- 
tional bounds are known, in the spirit of the work of H. Fujita and T. Kato 
[T9|J22j . to be sufficient for showing well-posedness for the 3D NS equations. 
Furthermore, they are directly linked to the helicity evolution for the dNS 
model, and therefore with a clear physical meaning and consequences. 
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1 Introduction 

The problem of whether the three-dimensional incompressible Navier-Stokes 
equations possess global regular and unique solutions is one of the most 
challenging unsolved problems in applied analysis [T51[Hi r24, 26 . In the two- 
dimensional case, global existence and uniqueness of weak and strong solu- 
tions is a well-established result; and in three dimensions, it is possible to 
show existence of weak solutions, globally in times, but we do not know how 
to establish their uniqueness for general initial conditions (see, e.g., [T5irHIl2~i] 
and references therein). Due to the notorious difficulties to attack the original 
problem in its whole generality, many sufficient conditions and criteria have 
been provided for establishing global regularity (for the most recent devel- 
opment in this regard see, e.g., [2l[5ll WTM2"Tll23U2{J] and references therein). 
Moreover, it may also be important and instructive to start from a decimated 
version of it, constraining the evolution of the Navier-Stokes velocity field in 
sub-manifolds with some given physical and dynamical properties. For ex- 
ample, it is well-known that viscous flows with certain physical symmetries, 
such as axi-symmetric flows [55] and helical flows [57J , have global regularity 
despite the fact that they have nontrivial 3D vortex stretching terms. Two 
new examples have recently appeared in the literature. In the first case |18j . 
the evolution of the Navier-Stokes equations is projected on a Fractal-Fourier 
space. In the second case [I], the velocity field is first exactly decomposed for 
each Fourier mode in two components, carrying positive and negative helic- 
ity respectively, and then evolved by keeping only one of the two, such as to 
make also the helicity invariant sign-definite (i.e. positive). Both approaches 
have been proposed in order to understand key phenomenological properties 
of turbulence as a function of its total number of degrees-of-freedom (for the 
Fractal-Fourier decimation) or as a function of the physical character of the 
non-linear interactions (for the helical-Fourier decimation; this is to distin- 
guish from the above mentioned flows with physical helical symmetry [161 
127] ). Clearly, it is tempting to ask also questions concerning regularities of 
their solutions. By investigating this model we hope to learn about possi- 
ble crucial physical features of the "long sought" singular solutions of the 
original Navier-Stokes equations that are eventually absent in our decimated 
version and in models with physical spatial symmetries. A similar approach 
has recently been proposed in |20j where the existence of global regularity for 
solutions of a modified Navier-Stokes equation has been proven. Let us notice 
however, that at difference from the previous cases, the equations studied in 
[20j cannot be seen as a decimation of the original Navier-Stokes case. 

In this paper we will concentrate on the regularity properties of the 
Navier-Stokes equations under helical-Fourier decimation only. Without loss 
of generality, we will study the case where the projection is made into the 
sub-manifold where helicity is positive-definite. It is worth mentioning that 
in this situation the positive-definite helicity is equivalent to the square of 
the H 1 / 2 — Sobolev norm. 

In such a case, we are able to show that the dNS system possesses global 
unique solutions which depend continuously on the initial data. Such a re- 
sult is obtained thanks to two additional estimates available only when the 
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dynamics is projected on a manifold with a sign-definite helicity. From a phe- 
nomenological physical point of view, this dNS equations also show a non- 
typical inverse-energy cascade, i.e. when forced on an intermediate range in 
Fourier space, energy tends to go to larger and larger scales and helicity is 
transported toward small-scales [I] . This is the so-called split-cascades regime 
typical of 2D Navier-Stokes flows - where enstrophy flux replaces the helicity 
(which is a positive-definite quantity here) flux. Such a phenomenology has 
also been observed in real 3D Navier-Stokes equations under fast rotation 
and with helical forcing [2"51l31| . 

The article is organized as follows. In section ([2| we review the main 
properties of the helical-decimated equations. In sec. ^ we briefly discuss 
the existence of weak solutions to the dNS system following the usual con- 
struction based on the Galerkin approximation procedure. In sees. ( 
establish two additional global uniform bounds on the velocity field, for initial 
data with finite positive helicity, by exploiting the constraints from the invis- 
cid conservation of helicity. Furthermore, we use these additional estimates 
to prove the uniqueness of the more regular solutions of the dNS system, i.e. 
solutions with initial data of finite positive helicity. Conclusions and a brief 
discussion on further perspectives can be found in sec. (rot. 



2 The helical-decimated Navier-Stokes equations 

The starting point of our analysis is the well-known helical-Fourier decom- 
position [32] of the velocity field v(x,t), expanded in Fourier series, u(k,t). 
Here we consider spatially periodic flows with zero spatial mean; consequently 
the zero wavenumber Fourier mode u(0,i) = 0. Being divergence-free, each 
velocity component in Fourier space has only two degrees of freedom. The 
idea is to define the two independent degrees of freedom by a projections on 
two orthonormal vectors, each one bringing a definite sign of helicity: 

u(k,t) = u+(k,t)h+(k) + u-(k,t)h-(k) (1) 

where h ± (fc) are the eigenvectors of the curl operator ikxh ± (k) — ±kh ± (k). 
In particular, one can always choose ft. ± (fc) = p,(k) x k ± ift, where fi is an 
arbitrary unit vector orthogonal to k which satisfies the relation p,(k) — 
—fi{—k) (necessary to ensure the reality of the velocity field). Such require- 
ment is satisfied, e.g., by the choice fi{k) — z x k/\\z x fe||, with z an arbi- 
trary vector. We then have for energy, E(t) = J d 3 x\v(x,t)\ 2 , and helicity, 
W(t) = J d 3 x v(x, t) ■ uj(x, t) (here uj(x : t) = V x v(x, t) is the vorticity): 

£(i) = E fc K(M)l 2 + K(M)l 2 ; m 
W) = E fc MK(M)l 2 -\u-{kM)- 1 1 

The non-linear term of the NS equations can be exactly decomposed in 4 
independent classes of triadic interactions, determined by the helical content 
of the complex amplitudes, u Sk (k) with Sk = ± (see [31] and Fig. [I]). Among 
three generic interacting modes u Sk (k) , u Sp (p) , u Sq (q) , one can identify 8 dif- 
ferent helical combinations (sk = ±, s p = ±,s 9 = ±). Among them, only 
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Fig. 1 Under helical decomposition the three wavenumbers interaction of the non- 
linear NS equa tion s are decomposed in four classes, depending on the relative he- 
licity signs. In 32 a simple dynamical argument is given supporting the fact that 
triads of classes (III) and (IV) mainly transfer energy toward small scales (high 
wave- numbers), i.e. they have the usual direct cascade, triads of class (I) enjoys 
an inverse energy cascade, while class (II) is mixed. In the figure this is summa- 
rized by red arrows denoting a backward energy transfer and by blue arrows for 
forward energy transfer. In [4] a direct numerical integration at high resolution of 
the dNS with only triads of class (I) showed that a stationary turbulent inverse 
energy cascade is indeed established. 



four are independent because of the symmetry that allows to change all signs 
of helicity simultaneously. 

Therefore, the Navier-Stokes equations in the helical- Fourier basis are [32] : 




fc+p+q=0 s p ,s q 

where v is the viscosity, f Sk (k), the external forcing and where 7 stands for 
complex conjugate. It is important to remark that the non-linear dynamics 
formally preserves both energy and helicity, triad-by-triad, i.e. dt(\u Sp \ 2 + 



\u s *\ 2 + \u Sk \ 2 ) = d t (s p p\u s "\ 2 + s q q\u s »\ 2 + s k k\u Sk \ 2 ) = 0. As a result, the 
evolution for helicity "formally" becomes: 



Among the 4 possible different choices of helicity signs interacting in the 
nonlinear term we will concentrate here on what happens when only one class 
is present. In particular, we will ask what happens when only interactions 
having the same sign of helicity for all the three wavenumber in each triad, 





j t U{t) = -^fc 3 (| U +(fc,t)| 2 - \u-(k,t)\ 2 ).+ 

k 

2Re^2k(u+{k,t)J + {k) - u~ (k,t)J~ (k)). 



(4) 



k 
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i.e. helicity becomes sign-definite. In order to do that we define the projector 
on positive/negative helicity states as 

v± HHk^hHk)^ 
h±(fc).fc±(fc) 

The projector is self-adjoint and commutes with derivatives. The action of 
this projector on a field in real space is defined via its Fourier decomposition, 
for example for projection on positive helicity we have: 

v+{x,t) = V + v{x,t) = ^e ik - x V + {k)u{k.t) =^e lk - x u+{k,t)h + {k). 

k k 

(6) 

We can then consider the dynamical evolution of the dNS: 

jd t v+ = V + {-v+ ■ Vv+ - Vp + ) + vAv+ + f+ 

where p is the pressure and / + is the external forcing, here taken for simplic- 
ity time independent and with zero mean such as we can restrict the analysis 
of the solutions of ([7]) in the space where J v + dx = 0. Following the tradition 
we denote by X the linear subspace of all functions in X that have spatial 
mean zero. Notice, that it is easy to realize that if the initial velocity configu- 
ration is chosen with only positive helicity components, i.e. v~(-,t — 0) = 0, 
the dynamical evolution of Q preserves this property for all times. 



3 Global existence of weak solutions with initial data in L 2 

The notion of weak solution for the dNS system ^ is similar to that of the 
3D NS equations [13,24 . Specifically, for a given initial data v + (-,t = 0) = 
Wp~(*) € L 2 , a weak solution of ^ in the interval [0, T] is a divergence-free 
vector field v+(x, t), which belongs to C weak ([0, T], L 2 ) f| L 2 ([0, T],H r ), and 
for which d t v+ e L 5 / 4 ([0, T], H' 1 ); such that the first equation of Q holds 
in I/ 5/,4 ([0, T], H^ 1 )- It is not difficult to prove the existence of weak solution 
for the dNS. It is enough to proceed as for the standard NS case, introducing 
a Galerkin projection on all wavenumbers smaller than a given value, N: 

v N (x,t)=V N v(x,t)= e lkx u{k,t), (8) 

|fc|<JV 

and then controlling the limit of the Galerkin approximates for N — > oo. 
The only additional detail consists in checking that the Galerkin projector, 
Vn commutes with the helical-Fourier projector V + , which is of course the 
case. We can then define a Galerkin-helical-Fourier projected field, evolving 
according to the finite system of ordinary differential equations: 



d t v+ = V+V N (-v+ ■ Vv+ - V P + ) + vAv+ + f+ 
V ■ v + - n 



(9) 

N = u - 
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Since the nonlinearity in the above ODE system is quadratic it has a short 
time existence and uniqueness solution for all initial data. As in the NS 
equations case one can show that the L 2 — norm of the solutions to the above 
system remain bounded and hence it has global existence and uniqueness. 
Following the same steps as for the standard original NS case [T5IrT4] starting 
from any finite N and then passing to the limit, modulo a subsequence Nj — > 
oo, using the relevant compactness theorem of Aubin |13| . 



4 More regular solutions with initial data in H 1 ^ 2 

In this section, we are going to use the fact that for the inviscid and unforced 
dNS ([9| system the helicity is formally conserved, and that it is positive- 
definite quadratic quantity, which is equivalent to the square of the H 1 / 2 — 
Sobolev norm. Therefore, obtaining uniform (in time) bounds on the helic- 
ity enables us to prove the existence of solutions with a higher degree of 
regularity, provided the initial data is in H 1 / 2 . Furthermore, this additional 
regularity will allow us to prove, in the next section, the uniqueness of these 
regular solutions within the class of weak solutions. Let us observe that all the 
estimates that follow are formal, but can be rigorously justified by obtaining 
them first for the corresponding solutions of the Galerkin approximating sys- 
tem ([9]), and then passing to the limit, modulo subsequences, with N — > 00. 
Furthermore, it is worth mentioning that similar ideas and estimates can be 
found in [19,22] in the study of short time existence and uniqueness of the 
three-dimensional NS equations with initial data in H 1 / 2 . The advantage of 
system ^ over the NS equations is that the H 1 / 2 remains finite, which al- 
lows to extend the short time existence argument to prove global regularity 
for all time and all initial data in H 1 / 2 . Indeed, the fact that helicity is a 
(positive-definite) inviscid invariant of |7]) can be readily verified, formally, 
by its evolution written in terms of the Fourier components Q , which after 
decimation reads: 

k k k 

where k = \k\, and u + (k,t) = T' + (k)u(k,t). Furthermore, the right-hand 



side (RHS) of ( 10 ) can be bounded by: 



RHS <J2\^u + (k,t)\\f+(k)k-i\ < E* 3 l« + (*.*)l 3 ) 4 El/ + (*)l 2 *" 1 ) 4 . 

(11) 

where the last step is obtained by using Cauchy-Schwarz inequality. By mul- 
tiplying and dividing by y/v and using Young's inequality we can further 
bound it by: 

(i/^fe 3 i« + (*,*)i 2 ) 4 Ei/ + ( fc )i 2 *~ 1,/ ~ 1 ) i ^ ( 12 ) 

k k 

^^(M^ + i-^l/Wfe- 1 . (13) 
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Plugging back the above estimate in (10) we get: 



1 d 
2dt 



£ k\u+(k, t)\ 2 +"-J2 * 3 l« + (*> t)\ 2 <lY, l/W*" 1 - (14) 



Observe that for every <f> g H s (i.e. belongs to the Sobolev space H s of 
periodic functions with zero average, J <f>(x) dx = 0) the H s — Sobolev norm, 
for s > 0, is equivalently defined by 



\\<t>\\Hs = (£k 2s \m\ 2 ) 1 



/2 



(15) 



Thanks to the above observation we can rewrite ( 14 ) in terms of the H 2 and 
i?s norms as: 



1 d 
2dt l 



v + \\ 2 1 +^> + ii 2 3 <iy;i/ + ( fc )i 2 * _1 - 



(16) 



Expression ( 16 1 can be rewritten in its equivalent integral form as: 

\\v + (t)\\ 2 Hh + f dt^\\v + (t')\\ 2 Hl < 11^(0)11^+ r^i-^iz+^oi^- 1 . 

(17) 

The above inequality implies that there is a weak solution of the dNS Q 
which satisfies: 



v+ e L°°([0,T];^); V + G L 2 ([0, T]; H*), 



(18) 



provided that / e £ 2 ([0, T]; £T - s) and v + (0) e iiX In the next section we 
will show that this more regular solution, that satisfies ( 18 ), is unique within 
the class of weak solutions. 



5 Uniqueness of more regular solution 

In order to prove the uniqueness of the more regular solutions with initial 
data in H 1 ' 2 , within the class of weak solutions, let us first rewrite the dNS 
equations Q in a symbolic form: 

d t v+ -i>Av+ = -B n (v+, v + ) + f+ (19) 

where B u {a,b) = P+[(aV) b - A~ l Wab)} 

Below, we give a formal proof of uniqueness which can be made fully rigorous 



by working with the integrated -in time- version of ( 20 ) and then following 
exactly the same argument as in Serrin [30] (see also [1|). It is worth men- 
tioning that similar ideas and estimates can be found in [H?I|2"2"] in the study 
of short time well-posedness of three-dimensional NS equations with H 1 / 2 
initial data. Let us begin by considering two different weak solutions, and 
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v^, with the same initial data in i/ 1 / 2 , and where satisfies (18). Their 
difference, w + = — v%, satisfies the evolution equation: 

d t w + - vAw + = -B-h(v+,w + ) - B H {w + ,v+) + B H (w + ,w + ). (20) 

Next, we study the time evolution of ||w + |||2 = / d 3 x\w + \ 2 , by formally 
multiplying (20) and integrating over the spatial domain. Note that V + is 
a self-adjoint operator, and that J d 3 x (V + a + ) ■ b + = j d 3 xa + ■ V + b + = 
j d 3 xa + ■ b + . By incompressibility, and after integration by parts and using 
the periodic boundary conditions, only the second term in the RHS of ( 20 ) 
survives and we formally obtain: 



w + \\l 2 +v\\Vw + \\l 2 =- J d 3 x[{w+ ■V)v+]-w+. (21) 



H 
2d~V 

The RHS of ( |2l| can be further integrated by part, using once again the 
incompressibility and the periodic boundary conditions, to get f d 3 x[(w + ■ 
V)i>i ] • w + = — J d 3 x[(w + ■ V)«; + ] • ; and by applying Holder's inequality 
one obtains the following bound: 

| J d 3 x{(w+ ■V)w+] ■ v+\ <\\w + \\ L 3\\Vw+\\ L2 \\vt\\ L ° (22) 

Applying the three-dimensional Sobolev embedding inequalities: 
II • ||l 3 < const. 1 1 • 1 1 #1/2 || • 1 1 lb < const. 1 1 • 



to (22 ) we obtain: 

| J d 3 x[(w+ ■V)w + ]-v+\<const.\\w + \\ H i\\Vw + \\ L 2\\v+\\ H i. (23) 

We now use the interpolation inequality || • < || • ||?j-«i|| ■ |lij*2 j with 
s = 9st + (1 - 9)s 2 with 6 = 1/2, Si = 0, s 2 = 1 to get: 

\\w + \\ Hh < Ih+HlolKllt - |k + ||| 2 ||Vu; + ||J 2 ; (24) 
and with 9 = 1/2, S\ = 1/2, s 2 = 3/2 to obtain: 

iKii^<iKiiiiiK + ni a - (25) 

n 2 M 2 



Putting all together we can bound the LHS of (22) as: 

^const.llV^+ll^ll^+ll^lKll^lKH^. (26) 
Let us now multiply and divide by v i : 



<const.[ V i ||V W +||| 2 ] B*.; (27) 
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and apply Young's inequality with the weights 3/4 and 1/4: 



const 



(28) 



Going back to equation (21 1 we can then bring the first term of the above 
relation on the LHS and get: 



1 d 1 ll™" 1 

\j t \\™ + \\h + \v\\Vw+\\l 2 <K — 



I L 2 I 



H 2 



(29) 



where we have used the assumption that € £°°([0, T]; ijs ), for every 
T > 0, since it satisfies (18). By virtue of Gronwall's inequality we finally 
arrive to: 



\w + (t)\\ 2 L2 < || W + (0)||i 2 exp(^- 3 ^ dt'\\vt(t')\\ 2 Hl ) (30) 



which proves uniqueness of the solution thanks to the fact that we have € 
L 2 ([Q,T}; H^) (since it satisfies (18 1. The above also implies the Lipschitz 
continuity with respect to the initia data due to the ||if + (0)||^2 prefactor. 



6 Conclusions and further perspectives 

Inviscid invariants of the NS equations are crucial in determining the physics 
of turbulence |17| . In 2D turbulence the presence of two positive-definite 
quadratic quantities, energy and enstrophy, leads to a split cascade: energy 
flows toward large scales and enstrophy to small scales p^lMTi?]. The fluid 
equations possess two inviscid invariants also in 3D: energy and helicity. At 
variance with energy, helicity is in general not sign-definite, and therefore it 
is not a 'coercitive' quantity. In principle, this allows for a simultaneous for- 
ward transfer of energy and helicity, as confirmed by the results of two-point 
closures [7,32 and direct numerical simulations p^0] J TT| . [28] . Nevertheless, a 
reversal of the flux of energy has been observed in geophysical flows subject 
to Earth's rotation |2"8"ll3"T] if the forcing is chiral, injecting a non-vanishing 
helicity in the system. For fast rotating turbulence in the inverse energy cas- 
cade regime, helicity is observed to still flow toward small scales: we have 
a split-cascade scenario as for the energy-enstrophy case of 2D turbulence. 
Further, the link between the existence of intermittent burst in the energy 
cascade and a "local" helicity blocking mechanism has been proposed [3]. 

In this paper we have shown that helicity plays also a very peculiar role in 
the determining the 'regularity' properties of the velocity field. If the dynam- 
ics is restricted to the sub-set of modes with a well definite sign of helicity (i.e. 
positive) then the flow admits unique global weak solutions that depend con- 
tinuously on the initial data. In our system the vortex stretching mechanism 
does not vanish, even if it is certainly depleted because of the local 'Bel- 
trami' condition imposed at each wave number by our decimation rule. Such 
a result leads to speculation that possible 'singular' solution of the original 
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NS equations should then be searched elsewhere, i.e. in the three classes of 
non-linear interactions connecting triads with different helicity components 
which arc killed in our decimated system. 
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